Search for diffuse neutrino flux from 
astrophysical sources with MACRO 



The MACRO Collaboration 

M. Ambrosio^^, R. Antolini^, G. Auriemma^'^'", D. Bakari^'^^, A. Baldini^^, 
G. C. Barbarino^^ B. C. Barisli^ G. Battistom6'^ Y. Bechermi^, R. Bellotti\ 
C. Bemporad^'^ , P. Bernardini^°, H. Bilokon^, C. Bloise^, C. Bower^, 
M. Brigida^, S. Bussino^^, F. Cafagna^, M. Calicchio^, D. Campana^^, 
M. CarboniS, R. Caruso^, S. Ceccliini^'^, F. Cei^^, V. Chiarella^, 
B. C. Choudhary^ S. Coutu"'% M. Cozzi^, G. De Cataldo^, H. Dekliissi^'i^, 
C. De Marzo\ 1. De Mitri^^, J. Derkaoui^-^^, M. De Vincenzi^^ 
A. Di Credico'', O. Erriquez^, C. Favuzzi^, C. Forti^, P. Fusco^, 
G. Giacomelli^, G. Giannini^^''^, N. Giglietto^, M. Giorgini^, M. Grassi^^, 
A. Grillo^, F. Guarino^^, C. Gustavino^, A. Habig^-P, K. Hanson^\ 
R. Heinz^, E. larocci^'*^, E. Katsavounidis^''^, I. Katsavounidis^''^, E. Kearns^, 
H. Kim^, S. Kyriazopoulou^, E. Lamanna^'^'', C. Lane^, D. S. Levin^^, 
P. Lipari^^, N. P. Longley^''*, M. J. Longo^^, F. Loparco^, F. Maaroufi^'^^, 
G. Mancarella^°, G. Mandrioli^, A. Margiotta^, A. Marini^ D. Martello^°, 
A. Marzari-Clliesa^^ M. N. Mazziotta^ D. G. Michael^, P. Monacelli^ 
T. Montaruli^ M. Monteno^^ S. Mufson^ J. Musser^ D. Nicol6^^ 
R. Nolty^ C. Orth^, G. Osteria^^^ O. Palamara^ V. Fateia^'", L. Patrizii^, 
R. Pazzi^^ C. W. Peck^ L. Perrone^"'* S. Petrera^ P. Pistillil^ V. Popa^'S, 

A. Raino^, J. Reynoldson^, F. Ronga^, A. Rrhioua^'^^, C. Satriano^^'*^, 
E. Scapparone'', K. Scholberg^''^, A. Sciubba^'^, P. Serra^, M. Sioli^, G. Sirri^, 
M. Sitta^^'°, P. Spinelli\ M. Spinetti^ M. Spurio^ R. Steinberg^ 
J. L. Stone^, L. R. Sulak^, A. Surdo^", G. Tarle", V. Togo^ M. Vakili^^'^ 

C. W. Walter^ and R. Webb^^ 



1. Dipartimento di Fisica dell'Universita di Bari and INFN, 70126 Bari, Italy 
2. Dipartimento di Fisica deH'Universita di Bologna and INFN, 40126 Bologna, 

Italy 

3. Physics Department, Boston University, Boston, MA 02215, USA 
4. California Institute of Technology, Pasadena, CA 91125, USA 

5. Department of Physics, Drexel University, Philadelphia, PA 19104, USA 

6. Laboratori Nazionali di Frascati dell'INFN, 00044 Frascati (Roma), Italy 

7. Laboratori Nazionali del Gran Sasso dell'INFN, 67010 Assergi (L'Aquila), Italy 
8. Depts. of Physics and of Astronomy, Indiana University, Bloomington, IN 

47405, USA 

9. Dipartimento di Fisica dell'Universita dell'Aquila and INFN, 67100 L'Aquila, 



Preprint submitted to Elsevier Science 



1 February 2008 



Italy 

10. Dipartimcnto di Fisica dell'Universita di Lecce and INFN, 73100 Lecce, Italy 

11. Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA 
12. Dipartimento di Fisica dell'Universita di Napoli and INFN, 80125 Napoli, Italy 

13. Dipartimento di Fisica dell'Universita di Pisa and INFN, 56010 Pisa, Italy 

14. Dipartimento di Fisica dell'Universita di Roma "La Sapienza" and INFN, 

00185 Roma, Italy 

15. Physics Department, Texas A&M University, College Station, TX 77843, USA 
16. Dipartimento di Fisica Spcrimentale dell'Universita di Torino and INFN, 

10125 Torino, Italy 

17. L.P.T.P, Faculty of Sciences, University Mohamed I, B.P. 524 Oujda, Morocco 
18. Dipartimento di Fisica dell'Universita di Roma Tre and INFN Sezione Roma 

Tre, 00146 Roma, Italy 
a Also Universita della Basilicata, 85100 Potcnza, Italy 
b Also INFN Milano, 20133 Milano, Italy 
c Also Istituto lASF/CNR, 40129 Bologna, Italy 
d Also Universita di Trieste and INFN, 34100 Trieste, Italy 
e Also Dipartimento di Encrgctica, Universita di Roma, 00185 Roma, Italy 
g Also Institute for Space Sciences, 76900 Bucharest, Romania 
h Macalester College, Dept. of Physics and Astr., St. Paul, MN 55105 
i Also Department of Physics, Pennsylvania State University, University Park, PA 

16801, USA 

I Also Dipartimento di Fisica dell'Universita della Calabria, Rende (Cosenza), 

Italy 

o Also Dipartimento di Scicnzc c Tccnologie Avanzate, Universita del Piemonte 

Oricntale, Alessandria, Italy 
p Also Duluth Physics Department, University of Minnesota, Duluth, MN 55812 
q Also Dept. of Physics, MIT, Cambridge, MA 02139 
r Also Intervideo Inc., Torrance CA 90505 USA 
s Also Resonance Photonics, Markham, Ontario, Canada 

* Corresponding author. Email: lorenzo.perrone@le.infn.it 



2 



Abstract 

Many galactic and extragalactic astrophysical sources are currently considered promis- 
ing candidates as high energy neutrino emitters. Astrophysical neutrinos can be 
detected as upward-going muons produced in charged-current interactions with the 
medium surrounding the detector. The expected neutrino fluxes from various mod- 
els start to dominate on the atmospheric neutrino background at neutrino energies 
above some tens of TeV. We present the results of a search for an excess of high 
energy upward-going muons among the sample of data collected by MACRO dur- 
ing ~ 5.8 years of effective running time. No signiflcant evidence for this signal 
was found. As a consequence, an upper limit on the flux of upward-going muons 
from high-energy neutrinos was set at the level of 1.7 xlO~^^cm~^ s~^ sr~-^. The 
corresponding upper limit for the diffuse neutrino flux was evaluated assuming a 
neutrino power law spectrum. Our result was compared with theoretical predictions 
and upper limits from other experiments. 

Key words: Neutrino, AGNs, GRBs 
PACS: 98.70.R 



1 Neutrino astronomy: overview and motivation 

Neutrinos with energies larger than ~ 1 GeV are expected from a wide class 
of galactic and extragalactic astrophysical sources. Neutrino production re- 
quires the existence of hadronic processes and it is generally described in the 
picture of the beam dum,p model [1]: high energy protons accelerated close to 
compact objects by shock waves or plasma turbulence interact with photons 
or target matter surrounding the source, producing pions. Neutrinos of elec- 
tron and muon flavors originate from decays of charged pions, as well from 
subsequent muon decays. In the same hadronic chains, high energy 7-rays are 
expected to be produced through neutral pion decays. Like 7-rays, neutrinos 
can travel undeficctcd through the Universe. Neutrinos however are much less 
absorbed than photons and thus make a more powerful probe for astronomy 
searches. Many of the candidate sources of neutrinos (binary systems, super- 
novae remnants, AGNs, GRBs, etc) have already been recognized as gamma 
rays emitters at energies higher than 1 TeV: this provides an important hint 
for neutrino astronomy, even though the observed 7-ray energies are not high 
enough to exclude the electromagnetic production mechanisms, such as syn- 
chrotron or inverse Compton processes. In this scenario, the detection of high 
energy neutrinos would open a new field of research, complementary to 7-ray 
astronomy. 

Neutrinos coming from a source propagate through the Earth, occasionally 
producing upward-going muons by charged-current interactions with the mat- 
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ter surrounding the detector. The detection probabihty grows with the energy 
due to the increase of both neutrino-nucleon cross sections and muon range, 
so that the effective detector mass results from the convolution of the detector 
area with the muon path-length. 

Underwater/ice neutrino detectors like Baikal [2] and AMANDA [3], as well 
as underground experiments like MACRO [4], 1MB [5], Baksan [6] and Su- 
perkamiokande [7], surveyed the sky in order to search for point source of 
astrophysical neutrinos. At the present their effective area is not large enough 
to measure the expected flux but they set experimental upper limits that 
guide and constrain theoretical models. Next generation neutrino telescopes 
like ANTARES [8], Icecube [9], NEMO [10] and NESTOR [11], with effec- 
tive areas ranging from 0.1 to 1 Km^, are expected to provide the sensitivity 
required for observing astrophysical sources. A detailed review dealing with 
physics goals and detection techniques in the held of neutrino astronomy can 
be found in [12]. 

The results of a search for point-hke sources was recently presented by the 
MACRO collaboration in [4] . Here we show the results of a search for a diffuse 
neutrino flux from unresolved sources, by looking for an excess of high energy 
upward-going events. A preliminary study of the sensitivity of the MACRO 
detector to high-energy muons is in [13]. In sec. 2 we discuss the capability 
of the MACRO detector as a neutrino telescope. In particular, since most of 
events induced by astrophysical neutrinos are expected in the 0.1-100 TeV 
energy range, the response of scintillation counters was accurately studied for 
high energy events, using the scintillator calibrations. In sec 3 we show that 
the energy released in the scintillation counters by crossing muons is the most 
important information for the selection of events induced by astrophysical neu- 
trinos. We describe the analysis cuts on the events from a simulated signal, 
and in sec. 4 on the events from the expected background. The main back- 
ground source are upward-going muons induced from atmospheric neutrinos. 
The GEANT-based [14] simulation tool used by MACRO [15] was modified 
to correctly handling the propagation of very high energy muons in the rock 
surrounding the detector [16], and the detector response for events with high 
energy releases. In sec. 5 we present the results of the analysis of the real data 
sample, and we give the upper limit for the diffuse neutrino flux from unre- 
solved sources. This upper limit is compared with some theoretical predictions 
and with upper limits from other experiments. 



2 MACRO as a neutrino telescope 

MACRO, located at the Gran Sasso Laboratories (Italy), was a large area 
multipurpose underground detector, in the shape of a rectangular box whose 
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global dimensions were 76.6 x 12.0 x 9.3 [15]. The lower half of the appa- 
ratus was filled with rock absorber alternating with streamer tube planes for 
particle tracking. Liquid scintillator layers, placed at the bottom, the center, 
the top and all around the detector provided time information for discriminat- 
ing the direction of incoming particles. The horizontal streamer tube planes 
were equipped with pick-up strips providing stereo read-out of the detector 
hits. The large acceptance (~ lO'* m^ sr for an isotropic flux), the good shield- 
ing of the site (the rate of cosmic ray muons is ~ 10~^ times the surface rate), 
the fast timing (time resolution ~ 0.5 ns) and the good pointing capabihty 
(intrinsic angular resolution < 1") made MACRO suitable for working as a 
neutrino telescope. 



2. 1 Scintillator response and photomultiplier saturation 

The energy deposited in a MACRO scintillation counter by a throughgoing 
particle is reconstructed by the ERP (Energy Reconstruction Processor) sys- 
tem according to the procedure described in [15]. 

The ERP circuit uses a twelve-bit ADC that makes an accurate integration of 
the photomultiplier pulse. This is applied to both the raw signal and a lOx at- 
tenuated version; the latter was used when the non attenuated ADC saturates, 
i.e. for fast magnetic monopole searches [17] and in this analysis. The ERP pro- 
vides also time of flight from two TDCs readout. The ERP ADCs and TDCs 
were calibrated weekly with delayed LED signals and variable-intensity laser 
pulses. The ADCs energy calibrations assume a linear relationship between 
the number of attenuated ADC counts and the light hitting the phototube L: 

L^G{ADC-p) (1) 

where G and p are the gain factor and the pedestal of each scintillation 
counter, respectively. This hypothesis is valid for most of the analyses based 
on MACRO data. 

The present analysis deals with signals greatly exceeding the average energy 
released in the detector by downward-going cosmic ray muons and atmospheric 
neutrino-induced events. If the energy released in a counter exceeds 4-=-5 times 
the level of that released by minimum ionizing cosmic ray muons (~ 34 MeV 
for one vertical track crossing a horizontal scintillation counter), the phototube 
response is not linear anymore [17]. The effect of photomultiplier saturation 
can be taken into account by including in eq. 1 the next-to-leading order terms. 
By following the approach described in [17], eq. 1 becomes: 

L = G{ADC -p) + Q{ADG - pf + T{ADG - pf + J={ADG - pY (2) 
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the constants Q, T, represent the second-, the third-, and the fourth-order 
corrections. These parameters were determined (one set for each MACRO 
scintillation counter) by fitting the laser light yield L versus the ERP [ADC — 
p) value with a fourth degree polynomial. G and p are fixed at the value given 
by the standard ERP energy calibration. 

In Fig. 1 (upper panel), the linear and the non-linear response regions are 
shown for one of the horizontal counters. The empty circles represent the mean 
values of the gaussian distributed ADC response for a fixed light level. The 
dashed curve shows the fourth-order fit done by taking into account the sta- 
tistical errors. The standard energy reconstruction is reliable below few times 
the MIP (Mnimum /onizing Particle) level. This analysis searches for events 
with a large energy release, above a certain threshold. For this reason, we used 
the correction given in eq. 2, which makes the energy reconstruction reliable 
above 3 MIP and up to ~ 15 MIP (i.e. ~ 500 MeV in a single counter). Above 
this value, the energy reconstruction may be affected by a large error because 
of the saturation. The fit constants Q, T, and the parameters G and p of 
eq. 2 were supposed independent on the time. This approximation is verified 
within 10% accuracy for the constants G and p measured with the standard 
calibration procedures. For the purpose of our analysis, a reconstructed energy 
> 500 MeV will be considered as an indication of a large energy release. 

The lower panel of Fig. 1 shows the distribution of the total energy released 
in the scintillators by a sample of 80000 real events collected by MACRO 
(downward-going cosmic ray muons). In particular, the solid line is obtained 
using eq. 2, whereas the dashed line is obtained using the standard calibra- 
tions (eq. 1). The corrections act only at large values of released energy Eji^i, 
while around the distribution maximum their effect is negligible. Ignoring the 
photomultiplier non-linearity leads to a systematic underestimate of the re- 
constructed energy loss for a small fraction of high energy events. This is in 
agreement with what it is expected from the curve given on the upper panel 
of Fig. 1 in the saturation region: a fixed ADC value corresponds to a value 
of the reconstructed energy higher than in the case of the linear response 
approximation. 



3 Monte Carlo simulation of the AGN signal 

This analysis aims to select a sample of very high energy upward-going muons, 
since at neutrino energies above some tens of TeV the predicted neutrino 

fluxes from astrophysical sources start to rise above the atmospheric neutrino 
background. In our simulation of the AGN signal in MACRO we used as input 
the model by Stecker et al. [18] for the photo-pion production of neutrinos in 
AGN core. 
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Fig. 1. (Upper) ERP attenuated- ADC response to laser light with linear and 
fourth-order fits. The plot is for one of the horizontal counters. The single {1 MIP) 
and ten-times {10 MIP) minimum ionizing particle levels are shown (for a vertical 
track crossing the center of the counter). (Lower) Distribution of the total released 
energy Er^i for a sample of 80000 downward-going muons in MACRO; dashed line: 
E^ei calculated with the linear response eq. 1; solid line: Eji^i calculated with eq. 2. 

Because high energy neutrino absorption from the Earth is an important effect, 
the neutrino propagation through the Earth was done by solving the kinetic 
equation for the transport of neutrinos in dense media. This approach is sug- 
gested in [19], where the neutrino absorption and the contribution of neutrino- 
nucleon neutral-current interactions were properly taken into account. The 
deep inelastic cross sections were calculated using the set of parton density 
functions CTEQ3-D1S [20]. The new version of parton density functions by 
CTEQ group [21] does not produce significant changes for our calculation. 
The neutrino-induced muon propagation in the rock outside the detector was 
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evaluated using the analytical formulas given in [22,16]. 

Then, following the analytical distributions, we simulated a sample of 13305 
upward-going muons on the surface of a volume containing the detector plus 
2 meter-wide layer of surrounding rock. The rock was included to take into 
account the effect of electromagnetic showers induced by the muons in the 
detector. In this volume, the transport of muons was done with the tool de- 
veloped in [16], which correctly propagate very high energy muons up to the 
PeV energies. This software package replaced in our simulation the default 
transport modules implemented in the GEANT [14] package. 

By normalizing to the expected event rate from the model by Stecker et al. 
(4.45 events/year on the surface of this box) our simulated sample is equivalent 
to Teg = 2988.5 years. 7547/13305 simulated muons reach the active detector 
and hit at least one scintillation counter in two different layers of the MACRO 
detector. This was the minimum requirement to define an "event". 

We adopt the following notation: taking as a reference the upper counter 
which measures the time Ti , the time of flight AT" — T2—T1 is positive if the 
particle travels downward and it is negative if the particle travels upwards. 
Two or more adjacent scintillator hits within a time window of 2 ns and on the 

same detector layer, at a maximum distance of 1 m, form a scintillator duster. 
Any association between a cluster and a single hit placed on different layers, as 
well as between two different clusters on different layers defines a scintillator 
track. The scintillator track length is the geometrical distance between the 
positions of the center of each cluster. If a cluster contains more than one 
counter, its center is calculated by averaging the hit positions weighted with 
the released energy. For each scintillator track, the quantity = cAT/L {L 
is the track length and c the speed of light) in our convention is around +1 
for downward-going particles and -1 in the opposite case. 

Fig. 2 (left) shows the correlation between the reconstructed total energy 
released in the scintillator Er^i, with the muon initial energy E initial- By initial 
energy we mean the energy of the muons when they enter the detector. Each 
simulated event is represented with a point in the scatter plot. 

In order to reduce the background due to atmospheric neutrino-induced muons 
(see next section), we select the Very High Energy (VHE) muons from the 
AGN simulated signal by imposing more and more stringent conditions, the 
effects of which are illustrated in Fig. 2 (right). The hght shaded histogram 
corresponds to events for which there is at least one scintillator with energy 
release greater than 500 MeV; the histogram with the dark shaded area corre- 
sponds to events with at least two scintillators, each one with energy release 
greater than 500 MeV, and distance between the hits smaller than 1 m; the 
histogram with the black shaded area corresponds to events with a further 
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Fig. 2. (Left) Scatter plot of the reconstructed total energy Ejiei released in the 
scintillator counters vs. the muon initial energy EjnUiai- Each simulated event is 
represented with a point. The darker points represent the subsample of events which 
survived all the cuts described in the text. (Right) The effect of more and more 
stringent cuts (see text) on Ej^i to select the subsample of high energy muons. 

scintillator in a different detector layer and with energy release greater than 
500 MeV. The darker points on Fig. 2 (left) represent the subsample of events, 
which survived all the cuts described above. Hereafter we will refer to the last 
condition as Ecut- In Fig- 2 (right) the principal maximum of the light shaded 
distribution is due to the muons which lose energy through standard ioniza- 
tion. The second "bump" is populated by higher energy muons, with large 
energy release in the scintillators through radiative processes. The value of 
ERei= 500 MeV was chosen after the calibration procedure described in sec. 
2. 



To discriminate the direction of the incoming events an algorithm based on 
the time-of-fiight technique was developed. The algorithm calculates the mean 
time associated with each scintillator layer involved in the event. We required 
a reconstructed upward-going flight direction between at least two different 
scintillator layers. This condition for selecting upward-going events was called 

D^cut- 

Many scintillator tracks (associations between any two scintillator clusters) 
are present in a showering event. Thus, a more selective cut for discriminating 
the upward direction, can be imposed using all the scintillator tracks in the 
event. We required that the ratio between the number of tracks with positive 
1/(3 value (downward-going) and negative 1/(3 (upward-going) is less than 0.4 
(cut D2cut)- The value of the ratio was optimized using a Monte Carlo study 
on a simulated sample of ~ 10^ downward-going muons: no downward-going 
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muons survived the energy cut Edt combined with the direction cut Dcut, 
defined as Dl^ut + D2cut- After the cuts Ecut + D^ut, 438/13305 simulated 
VHE neutrino-induced muons (signal) survived. 

Fig. 3 shows the distribution of the simulated cosine of muon zenith angle 
6sim- The sohd histogram refers to the overall sample of simulated events; the 
dashed histogram represents the sample selected after the Ecut + Dcut- The 
cos 9sim dependence of the muon flux is a consequence of high-energy neutrino 
absorption of the Earth [23]. Even assuming an isotropic neutrino flux [18] as 
input model, the zenith angular distribution of muons reaching the detector 
shows a strong suppression at cosOsim < —0.8 due to the hard discontinuity 
of the Earth density proflle in correspondence of the massive nucleus. 




-1 -0.8 -0.6 -0.4 -0.2 

Fig. 3. Distribution of the cosine of zenith angle for the overall sample of simulated 

events (solid line) . The dashed histogram refers to the events selected by the energy 
cut Ecut combined with the direction cut Dcut- 



Fig. 4 shows the selection efficiency of the analysis cuts Ecut + Dcut, when 
applied to the simulated sample of events from AGNs. In the left plot such 
efficiency is shown as a function of the cosine of muon zenith angle 9zen, while 
the right one represents the dependence on the muon initial energy. 

The energy cut Ecut and the flight direction cut Dcut were not completely effi- 
cient to reject atmospheric multiple muons and nearly horizontal high-energy 
cosmic ray muons. These events generate large showers in the apparatus, with 
a large number of scintillator tracks. Thus, a final cut was defined in order 
to improve the capability of selecting upward-going events against the back- 
ground of downward-going showering events. For each event, we define as: itot, 
the number of scintillator tracks whose length is greater than 2.5 m; iup, the 
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Fig. 4. Efficiency of analysis cuts E^ut + D^ut for selecting muons from the simulated 
AGN signal, as a function of the muon zenith angle (left) and of the muon initial 
energy (right). Only the statistical errors due to the limited simulation sample are 
shown. 

number of scintillator tracks for which —1.25 < 1//3 < —0.75; idown, the num- 
ber of scintillator tracks for which 0.75 < 1//3 < 1.25. The final cut (called 
Scut' itot > 11 and iup/itot > 0.3 and idown/iup < 0.1) selects upward-going 
events for events with a large number of scintillator tracks. This cut was 
tuned to reject Monte Carlo simulated background events and the large sam- 
ple of multiple muons and nearly horizontal high-energy cosmic ray muons 
in the real data. After this last cut, 247/13305 simulated VHE signal events 
survived. 



4 The background of atmospheric neutrinos 



The background from atmospheric neutrinos was estimated with the Monte 
Carlo program developed for the study of the atmospheric neutrino-induced 
upward-going muon flux in MACRO. It uses the Bartol neutrino flux [24], 
the Morfin and Tung parton density functions [25] for the calculation of the 
DIS uN cross sections, and the Lohmann et al. [22] muon energy loss for the 
propagation of induced muons through matter. The theoretical uncertainty 
on the expected event rate in MACRO for atmospheric neutrinos is ~ 17%, 
mainly due to the uncertainty on the neutrino flux. Further details can be 
found in [26-28]. 

Fig. 5 (left) shows the distribution of the muon initial energy (i.e. muon energy 
when entering the detector) for the atmospheric neutrino simulated sample 
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(dashed line) and for the AGN sample (solid line). The events surviving the 
analysis cuts E,.ijt + D^ut are also shown (light shaded histogram for the atmo- 
spheric sample, dark shaded histogram for the AGN sample). The Fig. 5 (right) 
shows the distribution of the reconstructed energy released in the scintillator 
counters, with the same notation of the previous plot. The two distributions 
(atmospheric neutrino background and AGN signal) have been normalized to 
the same equivalent time of the AGN sample. A small number of background 
events survive the analysis cuts Ecut + D^ut- Because of the Monte Carlo gen- 
eration statistics, the error bars (not shown) of the light shaded histograms, 
are at the level of 50% for each energy bin. 

As a final remark, it should be noticed that the simulation of the atmospheric 
neutrino background does not take into account the isotropic contribution 
of neutrinos from semileptonic decays of charmed hadrons, usually called 
"prompt neutrinos". For neutrino energies above 10 TeV [30-32] this con- 
tribution should begin to dominate over the conventional pion and kaon decay 
induced neutrino fiux. Actually, mainly due to the lack of precise informa- 
tion about high energy charm inclusive cross sections, the expected fluxes of 
prompt neutrinos, as well as the energy above which they start to dominate, 
fluctuate on a very wide range. Even by assuming one of the highest prompt 
neutrino flux prediction given in [32], the contribution to the signal has been 
estimated as few percent of the conventional atmospheric neutrino background 
(at least in the energy range in which this analysis is sensitive). 

In Fig. 7 the contribution from prompt neutrinos, according to the maxi- 
mum (PNmaa;) and the minimum (PN^in) predictions given in [32], are shown, 
together with the flux of conventional atmospheric neutrino and the diffuse 
fluxes of neutrinos from AGNs and GRBs according to some predictions. 



5 Data analysis and results 

The data used for this analysis were collected in the period from April 1994 to 
December 2000 (5.8 years, including efficiencies). After imposing the energy 
cut, Ecut-i and the flight direction cut, Dcuti 97 real events survive. 
Three events were rejected due to the action of a cut on the quality of the 
data. This quality cut (applied in all MACRO analysis which uses the ERP 
TDCs) requires that the two TDCs, which have a different dynamic range, 
on each ERP channel must agree. This inefficiency, calculated from different 
analysis, affect 2% of the data sample [27,28]. 

As expected, a large contamination from multiple muons and near horizontal 
events is found. Those background events were rejected with the cut Scuu after 
applying this cut, only two events survive. Fig. 6 shows the distribution of the 
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Fig. 5. (Left): distribution of the muon initial energy for the AGN sample (solid line) 
and for the atmospheric neutrino background (dashed line). The events surviving 
Ecut + Dcut are shown as the light shaded histogram (atmospheric sample) and as 
the dark shaded histogram (AGN sample). (Right): distribution of the reconstructed 
energy released in the scintillators, with exactly the same notation as for the plot 
in the left. 



quantity 1/P for one of the two selected event. The for the scintillator 
tracks fulfilling the conditions on geometrical length and on timing quality 
are shown as the shaded histogram. 

An attempt to reconstruct the direction of both events has been performed 
by combining scintillator counters with streamer tubes tracking information. 
Because of the strong showering activity, the result is affected by a too high 
uncertainty for any directional astronomy study. 

The simulated AGN signal events have been processed with the same analysis 
chain as the real data; they are reduced by a factor (2%) which takes into ac- 
count the above ERP inefficiency, not explicitly reproduced in the simulation 
of the detector response. 242/13305 simulated upward-going muons from dif- 
fuse AGNs survived all cuts. These events correspond to 0.54 for 5.8 years live 
time. For the atmospheric neutrino background, 1.1 events are expected in the 
same running period. The results are summarized in Table 1. The MACRO 
display of the two real events is shown in Fig. 8. 
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Fig. 6. 1/(5 distribution for one of the two events (Run 16399, Evt 2925) surviv- 
ing all the analysis cuts. One entry corresponds to one scintillator track. The 
scintillator tracks fulfilling both the conditions on geometrical length and on tim- 
ing quality are shown as the shaded histogram. 

5.1 The upper limit 



The data show no evidence of possible excess due to a diffuse neutrino flux 
from AGNs within the framework of the model [18]. Our result was used to 
set a muon flux upper limit Fj^ which was calculated at 90% confldence level: 



pL _ Upper limit (90%c.L) 

The calculation of the numerator is based on the observation of 2 events when 
the expected background is 1.1 events, according to the recent approach of 
Feldman and Cousins [33]. e is the fraction of simulated AGN events which 
survive the analysis cuts; Area{VL) is the geometrical area seen by the diffuse 
flux as a function of the solid angle Vt (the integral extends to the lower Earth 
hemisphere) and Ti — 5.8 years is the considered MACRO live time. With 
these values the muon flux upper limit becomes: 



F^(90%c.L) = (1.7 ±0.2) X lO-^^cm-2 s"^ sr 



(4) 



This limit can be finally converted into a (differential) neutrino fiux upper 
limit. In order to compare our results with the upper limits given by other 
experiments, it is convenient to assume a power law spectrum with spectral 
index 2 for the initial neutrino flux. We have then weighted the number of 
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T?atp of Riirvivpfl pvpnt.R in Pi 8 v 


Atmospheric (MC) 
Agn (MC) 
DATA 


1.1 ±0.5 
0.54 ±0.03 
2 



Table 1: Expected rate of events surviving all analysis cuts in 5.8 years of MACRO 
running. Besides the number of events in the data, the expectation from the atmo- 
spheric background and the AGN neutrino flux (according to [18]) are shown. Only 
the statistical errors of the simulation were considered. 

surviving events with the ratio between this spectrum and the spectrum cal- 
culated by Stecker et al. With these hypotheses, the neutrino flux upper limit 
from this analysis is ■ ^ (4.1 ± 0.4) xlO~*^ GeVcm"^ s~^ sr~^ Actually, 
this value must be considered carefully because it is affected by model as- 
sumptions (initial spectrum, neutrino cross sections and muon energy losses). 
Here, it has been estimated for comparison purposes. 

Fig. 7 shows some theoretical predictions of neutrino fluxes from astrophysical 
sources (AGNs and GRBs) in comparison with upper limits obtained by cur- 
rent experiments. A list of them with references and energy range of sensitivity 
is given in Table 2. 





■ (90% C.L.) 

(G(-Vcm-2s-isr-i) 


Neutrino Energy Range 

(GcV) 


EAS-TOP [37] 


2.0 xlO-3 


10^ 10^ 


SPS-DUMAND [38] 


6.0 xlO-^ 


10^ ^ 10^ 


Baikal [2] 


1.4 xlQ-^ 


10*^ ± 10^ 


Baikal (ue) [39] 


1.3 ± 1.9 xl0~6 


10^ ± 10^ 


Prejus [40] 


5.0 xlO"^ 


~ 2.6x10^ 


MACRO (this analysis) 


4.1 ± 0.4 xlO-6 


10^ 10^ 


AMANDA [41] 


1.0 xlO-6 


10^ 10^ 



Table 2: Neutrino flux upper limits (90% C.L.) from current experiments. The neu- 
trino sensitivity range is shown in the last column. 



15 




Fig. 7. Diffuse fluxes of t'^ + from AGNs and GRBs according to many predic- 
tions. (SDSS) Stcckcr et al, 1991 [18]; (M) Mannheim, 1995 [34]; (W) Waxman and 
Bahcall, 1998 [35]; (D) De Paolis et al, 2001 [36]. The dashed curve refers to the an- 
gle average atmospheric neutrino flux; (G) Gaisser et aL,1995 [1]. The dash-dotted 
hues refer to the flux of prompt neutrinos according to the maximum (PN^ax) and 
the minimum (PN^j„) predictions given in Costa, 2001 [32]. The dotted hue is the 
theoretical upper bound to neutrino fluxes from astrophysical sources as calculated 
in [35]. Some current published experimental upper limits are also shown. See Tab. 
2 for references to the experiment results. 



6 Conclusions 



We presented the results of a search for a diffuse neutrino flux from unresolved 
astrophysical sources by analyzing a sample of high energy events collected by 
MACRO. We tested the detector response at high energy taking special care 
for the reliability of the scintillation counter response to high-energy particles. 
We used the reconstructed energy released in the scintillators to select a sample 
of high energy events. Two high energy upward-going candidates have been 
found in 5.8 y of effective running time. This number is compatible with the 
expected atmospheric neutrino background (1.1 events). The analyzed data 
have been used to set a muon and neutrino flux upper limit comparable with 
the results given by other experiments. 
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RUN 16399 EVENT 2925 at 1:28:12.77 on 9/14/1998 
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RUN 19574 EVENT 2198 at 10:33:47.90 on 5/24/2000 




Fig. 8. Display of the two surviving events (longitudinal view) along the x-axis. The 
horizontal lines represent the 10+4 planes of horizontal streamer tube wires in the 
bottom and in the "Attico" parts of the detector. The wire hits are represented 
by black points. The gray boxes represent scintillator hits. 32 scintillation counters 
overall fired in the first event which correspond to 381 scintillator tracks, whose 
1/P values are shown in Fig. 5. 12 counters fired in the second event. The location 
of fired scintillation counters is also given: T— >• Top, Central, B— > Bottom, E— 
East, W— > West, North, S— > South, respectively. 
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